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Abstract 

In a seminal paper [12] , Weitz gave a deterministic fully polynomial approximation scheme for count- 
ing exponentially weighted independent sets (equivalently, approximating the partition function of the 
^si hard-core model from statistical physics) on graphs of degree at most d, up to the critical activity for 

the uniqueness of the Gibbs measure on the infinite d-regular tree. More recently Sly [10] showed that 
this is optimal in the sense that if there is an FPRAS for the hard-core partition function on graphs 
of maximum degree d for activities larger than the critical activity on the infinite d-regular tree then 
NP = RP. In this paper, we extend Weitz's approach to derive a deterministic fully polynomial approx- 
imation scheme for the partition function of the anti-ferromagnetic Ising model with arbitrary field on 
graphs of maximum degree d, up to the corresponding critical point on the d-regular tree. The main 
ingredient of our result is a proof that for two-state anti-ferromagnetic spin systems on the d-regular tree, 
weak spatial mixing implies strong spatial mixing. This in turn uses a message-decay argument which 
extends a similar approach proposed recently for the hard-core model by Restrepo et al [5] to the case 
1 of the anti-ferromagnetic Ising model with arbitrary field. By a standard correspondence, these results 

£_) translate to arbitrary two-state anti-ferromagnetic spin systems with soft constraints. 

1 Introduction 

1.1 Background 

fvj Spin systems are a general framework for modeling nearest-neighbor interactions on graphs, and are widely 

studied in both statistical physics and applied probability. A spin system consists of a large collection of 
nodes, each of which may be in one of a fixed number of states called spins. A neighborhood structure 
is specified by edges between the nodes. Interactions between neighboring nodes are determined by edge 

t-H potentials, which assign an energy value to each edge based on the spin values of its endpoints. In addition, 

there are vertex potentials which assign an energy value to each node based on the value of its spin. For any 

• ^h configuration a of spins on the nodes, the energy H (a) is just the sum of its edge and vertex energies. Based 

on the Gibbs formalism from statistical physics, the probability of finding the system in configuration a is 
then proportional to the weight w(cr) = exp(—H(a)). 

In this paper, we concentrate on two-state spin systems, where each vertex can be in one of two states, 
referred to as "+" and "— " . Such a system can be defined by specifying a (+, +) edge activity (3, a (— , — ) edge 
activity 7, and a vertex activity A, where /3, 7 and A are non-negative parameters. For a graph G = (V, E), 
a configuration a : V *— > {+ , — } is an assignment of + and — spins to the vertices of 67. The weight w(a) 
of the configuration a is given by 

w(a) = A" 1 ^)^^™-^), (1) 
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where m(a) denotes the number of vertices assigned spin — , and n + (a) (respectively, n_(<x)) denotes the 
number of edges for which both endpoints are assigned spin + (respectively, — ). The partition function of 
the model is defined as 

Z= J2 «»(*)■ ( 2 ) 

^e{+,-} v ' 

The partition function, in addition to being a natural weighted generalization of the notion of counting, 
is a fundamental quantity in statistical physics. For example, it is the normalizing factor in the Gibbs 
distribution: the probability of occurrence of configuration a is given by /i(cr) — w(a)/Z . In addition, many 
other properties of the model can be deduced by studying the partition function [3] . 

As a simple concrete example of a two-state spin system, consider the setting f3 = 1 and 7 = (so that 
configurations with adjacent "— " spins are assigned weight zero, and thus prohibited), known as the hard- 
core model. The associated Gibbs distribution is a weighted measure on independent sets in the graph G, 
in which any independent set U has weight X' u '. Another important class of examples, known as the Ising 
modeFj is obtained by setting f3 ~ 7 > 0. There is a significant qualitative difference between the Ising 
model with f3 = 7 > 1 (the ferromagnetic case) and with f3 = 7 < 1 (the anti-ferromagnetic case). The 
latter is an example of a "repulsive" model, which means that the edge potentials assign higher weights to 
edges with different spins at their endpoints, while the ferromagnetic case is "attractive" (higher weights 
are assigned to edges with the same spin at their endpoints). The parameter A can be identified with an 
"external field", i.e., a bias associated with each spin. The case A = 1 corresponds to zero field, while 
A < 1 and A > I correspond to positive and negative fields respectively. More generally, we will refer to any 
two-state system satisfying (3"/ > 1 as ferromagnetic, and any satisfying ^7 < 1 as anti-ferromagnetic. Also, 
a model satisfying ,87 > is said to have soft constraints (in the sense that no combination of spin values at 
adjacent vertices is prohibited). In a sense to be made precise later (see Appendix |A|, Ising models capture 
arbitrary two-state spin systems with soft constraints; in particular, the two descriptions are equivalent on 
regular graphs. For this reason we will henceforth focus mainly on Ising models. 

The theory of spin systems derives in large part from considering the limiting behavior of the Gibbs 
distribution as the size of the underlying graph goes to infinity. Based on the above formalism for finite 
graphs, one may define a Gibbs measure /i on an infinite graph Q by requiring that the marginal distribution 
on any finite subgraph H, conditional on the configuration on G\H, is given by equation (fTl). (Here the spins 
in Q\H act as a fixed boundary condition in (fTl). ) It is a well known result in the statistical physics literature 
(see, for example, [3]) that at least one such measure fi can always be defined. However, for certain values 
of the parameters of the spin system there may be multiple solutions for [i, in which case the Gibbs measure 
is said to be non-unique. 

We will now look at the phenomenon of non-uniqueness more closely in the special case when the infinite 
graph Q is a d-ary tree|j As noted above, the anti-ferromagnetic Ising model captures all two-state anti- 
ferromagnetic spin systems with soft constraints on regular graphs, and hence it is sufficient to consider the 
Ising case. Consider an anti-ferromagnetic Ising model on the d-ary tree with edge activity /?(= 7) and vertex 
activity A. It turns out that if /3 > ^^ then the Gibbs measure is unique for all values of A. In particular, in 
the zero-field case A = 1, the Gibbs measure is unique if and only if (3 > jt^- However, when (3 < 4=1 the 
Gibbs measure is no longer unique for all values of the vertex activity A. In this case, there exists a critical 
activity \ c ((3,d) > 1 such that the Gibbs measure is unique if and only if |logA| > log X c (f3, d). We sketch 
the curves of log A c (/3, d) in Figure [l] for d — 5 and d = 13. The area below the curves is the non- uniqueness 
region. We note here that the non-uniqueness region is monotonically increasing with degree, so the curve 
for d — 5 lies strictly below that for d = 13. Also, note that the curves intersect the /3-axis at /3 = |xj- 



1 The description of the Ising model given here differs slightly from the more popular description in terms of edge and vertex 
potentials outlined in the first paragraph. However, translating between the two descriptions is easy; see Appendix |A| 

2 We remark here that the infinite (d+ l)-regular tree and the infinite d-ary tree show exactly the same behavior with respect 
to the uniqueness of the Gibbs measure. This follows immediately from the fact that the (d + l)-regular tree can be viewed 
as a root attached to the roots of d + 1 infinite d-ary trees. We shall thus move freely between these two objects for ease of 
exposition throughout the paper. 
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Figure 1: logA c (/3, d) for d = 5 and d = 13. The curves intersect the /3-axis at j3 = | and /? = f respectively. 



The phenomenon of non-uniqueness of the Gibbs measure can also be described in terms of the more 
algorithmic notion of decay of correlations. We stick to our example of the infinite d-aiy tree. Fix a vertex 
v in the tree, and let Si be the set of vertices in the tree at distance at least I from v. Let q v {l,a) be the 
probability of having spin 4- at v conditional on the configuration on Si being a. It turns out that uniqueness 
of the Gibbs measure is equivalent to the condition that the inequality 



\q v (l,cr) - q v (l,T)\ < exp(-fi(Z)) 



(3) 



holds for any two configurations a and r on Si [j The above condition is referred to in the literature as weak 
spatial mixing. 

It has been believed for a long time (and proved in various manifestations) that there is an intimate 
relationship between weak spatial mixing and the running time of algorithms for approximating the associated 
partition function: roughly speaking, in the uniqueness region (where there is decay of correlations), the 
system should be amenable to local algorithms and thus be computationally tractable. A spectacular result 
in this direction was Weitz's fully polynomial deterministic approximation scheme (FPTAS) for the partition 
function of the hard-core model, which works on all graphs of degree at most d + 1 for all activities A less 
than the critical activity X c (d) for the uniqueness of the Gibbs measure on the infinite <i-ary tree [12 . This 
is even more interesting in light of a recent breakthrough due to Sly [TO] (see also [2]), who showed that the 
existence of an FPRAS for the partition function of the hard-core model on graphs of degree at most d + 1 
for activities larger than X c (d) would imply that NP = RP. Thus the range of validity of Weitz's algorithm 
is optimal. 

Weitz [12} gave a general two-step framework for designing deterministic algorithms for approximating 
partition functions of two-state spin systems. To describe this framework, we begin with the standard 
observation that in order to get an FPTAS for the partition function, it is sufficient to give an FPTAS 
for the probability of having spin + at any given vertex v. The first component of the framework is a 
combinatorial reduction, which shows that the problem of approximating this probability for a general two- 
state spin system on a graph G of maximum degree d + 1 can be reduced to the problem of approximating 
the same probability on a related finite subtree of the infinite (d + l)-regular tree rooted at v, in which the 
spins of some of the vertices are fixed to certain values (this is the so-called self-avoiding walk tree of the 



3 To be precise, this condition does not hold on the boundary of the uniqueness region, that is, for | log A | = logA c (/3,d); at 
this critical value, the l.h.s of equation J31 still decays to with I, but not at an exponential rate. We will focus on the interior 
of this region, and by a slight abuse of terminology refer to it as the "uniqueness region" . 



graph G). We emphasize that this is a model-independent reduction, and depends only upon the fact that 
the number of spin values is two. The associated self-avoiding walk tree, however, may be exponential in 
the size of the original graph G, and thus one needs to show that it is sufficient to truncate the tree at a 
depth logarithmic in the size of G in order to obtain a good approximation. However, since some of the 
fixed vertices in the tree might be very close to the root v, it is not possible to argue using weak spatial 
mixing that a logarithmic depth of recursion suffices for approximating the partition function (because the 
parameter I in equation (13J) must be taken to be the minimum distance of a fixed vertex from the root). 

Accordingly, the second component of Weitz's framework is to establish that, for the spin system in 
question, weak spatial mixing on the infinite ci-ary tree is in fact equivalent to strong spatial mixing, which 
roughly states that the exponential decay of point-to-set correlations (pi) guaranteed by weak spatial mixing 
holds also when the spins at an arbitrary set of vertices are fixed to arbitrarily chosen values (see Section [2] for 
a precise definition). Weitz [T2] established this fact for the hard core model, using a step- by-step comparison 
of ratios of occupation probabilities on the standard ci-ary tree and on the modified tree with fixed vertices. 
It was claimed in [12] that such a result holds also for the anti-ferromagnetic Ising model, but to the best of 
our knowledge no proof of this fact (except in the special zero-field case where A = 1; see [13]) has so far 
been published. 

1.2 Contributions 

In this paper, we give a proof of the fact that for the anti-ferromagnetic Ising model with any field, weak 
spatial mixing implies strong spatial mixing on the ci-ary tree. Formally, we have the following theorem. 

Theorem 1.1 For the anti-ferromagnetic Ising model with arbitrary field on the d-ary tree with d>2, weak 
spatial mixing implies strong spatial mixing. 

Notice that it is easy to see that this holds also for the infinite (d+ l)-rcgular tree, since the (d+ l)-regular 
tree and the ci-ary tree differ only in the degree of the root. We also note that by the translation described 
in Appendix [A] Theorem |1.1| holds also for arbitrary anti-ferromagnetic two-state spin systems with soft 
constraints. 

Given Weitz's general reduction described above, we obtain as an almost immediate consequence of The- 
orem |1.1| an FPTAS for the partition function of the anti- ferromagnetic Ising model on graphs of maximum 
degree at most d + 1 throughout the uniqueness region of the Gibbs measure on the d-ary tree. 

Corollary 1.2 Let d > 2. Consider an anti-ferromagnetic Ising model with parameters f3 and A. For /3 
and A in the interior of the uniqueness region of the d-ary tree, there is a deterministic polynomial time 
approximation scheme for the partition function of the associated spin system on graphs of degree at most 
d+l. 

By the translation described in Appendix[S] we can extend this result to general two-state anti-ferromagnetic 
spin systems. The difference is that the critical activity may now differ for vertices of different degrees. Let 
A c (/3, d) be the critical activity for the anti-ferromagnetic Ising model described above (and defined formally 
in Section|2.2.1). Then, we have the following corollary. 



Corollary 1.3 Let d > 2. Consider an anti-ferromagnetic two-state spin system with parameters /3, 7 and 
A. Let f3' be the edge potential for the equivalent anti-ferromagnetic Ising model. Let Q be the class of graphs 
with maximum degree d + 1 in which every vertex v satisfies the condition 

I log A, I 4 I log A + ^|^ (log 7 - log (3) I > log A c (/3', d). 

Then there is a deterministic polynomial time approximation scheme for the partition function of the associ- 
ated spin system on graphs in the class Q . In particular, the class Q includes all (d+ \)-regular graphs when 
/?, 7 and A are in the interior of the uniqueness region of the d-ary tree. 



We briefly sketch the approach we use to prove our main technical result, namely that weak spatial 
mixing implies strong spatial mixing (Theorem |1.1[ ) . Inspired by recent work of Restrepo et al [§] , we 
design a "message" (i.e., an invertible function of the probability of a vertex having spin + ) such that 
"disagreements" in the message decay by a constant factor at each vertex of the tree. The challenge is to 
ensure that such a message can be designed for all points in the uniqueness region of the d-axy tree. For the 
special zero-field case of the anti-ferromagnetic Ising model (when A = 1), such a message is well known 13 . 
However, this message does not work up to the threshold for general vertex potentials A. Restrepo et al [5] 
recently derived a message which works up to the tree uniqueness threshold for the hard-core model. For 
the general anti-ferromagnetic Ising model, such a message turns out to be more complex than those known 
for the zero-field case and for the hard-core model. Our message is defined at the beginning of Section [3] 
and the requisite decay property is established in Section [4] 

We conjecture that our proof of strong spatial mixing based on stepwise decay of messages may lead to 
further consequences. For example, as shown by Restrepo et al [9), the message decay property can be used 
to extend Weitz's algorithm by exploiting the structure of special classes of graphs to obtain approximation 
algorithms beyond the tree threshold for those graphs. In addition, our proof demonstrates the versatility 
of the message approach. 

Finally, we point out a byproduct of our approach that may be of independent interest. The exact 
value of the critical field A c (/3, d) as a function of (/3, d) is apparently widely accepted folklore, but the only 
derivation we could find for it in the literature [3J p. 255] does not provide a formal proof, appealing instead 
to numerical evidence. We point out that our proof of strong spatial mixing does not assume knowledge 



of X c (/3,d), and in fact derives its value as a byproduct (see the Remark following Theorem 2.5 for more 
details). Thus our approach gives a proof of the location of the uniqueness threshold A c (/3, d) of two-state 
anti-ferromagnetic spin systems. 

Remark: After obtaining our message-decay proof, we received a sketch of Weitz's original unpublished 
proof [11 . It is interesting to note that that proof is quite different from ours, and employs a delicate 
two-step analysis of the tree recursion described in Section [2J For reasons mentioned above, we believe that 
our message-decay proof, in addition to being the first published version of this result, is potentially more 
robust and flexible than Weitz's approach; for example, it is not clear how to adapt Weitz's analysis to 
obtain stronger results for special classes of graphs such as lattices, as is done in [S]. 

1.3 Related work 

Our work is mainly motivated by the deterministic counting algorithm of Weitz |12j , which was the first to 
show an interesting connection between the running time of an algorithm not related to Markov chain Monte 
Carlo and the phase transition phenomenon for spin systems. On the complexity side, using a randomized 
gadget first proposed by Dyer, Frieze and Jerrum pQ and analysed further by Mossel, Weitz and Wormald [8], 
Sly jTD] proved that if there is an FPRAS for the partition function of the hard-core model on graphs of 
degree at most d in the non-uniqueness region of the d-regular tree, then NP = RP, thus showing that the 
range of validity of Weitz's algorithm is optimal. Technically Sly's result holds only sufficiently close to the 
boundary of the uniqueness region; this restriction was mostly removed in a recent paper of Galanis et al [2] . 
For the case of unbounded degree graphs, Goldberg, Jerrum and Paterson [5] showed that approximating 
the partition function for the zero-field case (A = 1) is NP-hard in the interior of the square < /3,7 < 1. 

A related problem is to get exponential lower bounds on the mixing time of any local Markov chain 
(Glauber dynamics) that samples from the hard-core and anti-ferromagnetic Ising models. Mossel, Weitz 
and Wormald [8J and Gerschenfeld and Montanari |3] showed that beyond the uniqueness threshold for 
d-regular trees, Glauber dynamics for these models can take exponential time to mix on d-regular graphs. 
Gerschenfeld and Montanari also show that for these models on random regular graphs, this threshold for 
slow mixing is also the threshold beyond which the reconstruction problem is solvable, and pointed out that 
these results therefore establish the existence of (i-regular graphs on which the reconstruction problem is 
solvable beyond the uniqueness threshold for the (i-regular tree [I]. 

On the algorithmic side, an analysis of Weitz's algorithm for the zero-field case of the anti-ferromagnetic 



Ising model appears in [T3]. There has been some subsequent progress on the hard-core model on special 
classes of graphs too: recently, Restrepo et al \§\ used a message-decay proof to get improved strong spatial 
mixing thresholds on the 2D integer lattice for the hard core model. They achieved this by exploiting the 
special structure of self-avoiding walk trees obtained when Weitz's reduction is applied to the lattice. The 
message-decay proof turns out to be crucial in tightening the analysis to obtain strong spatial mixing over a 
wider range of parameters for these special trees. Much more is known about algorithms for the ferromagnetic 
case: Jerrum and Sinclair [5] gave an FPRAS for the Ising model with arbitrary field on graphs of arbitrary 
degree, while Goldberg, Jerrum and Paterson [5] showed how to extend this to the whole of the ferromagnetic 
region (3j > 1 with A = 1. The latter paper [5] also gave an FPRAS for the partition function on graphs of 
arbitrary degree for parts of the anti- ferromagnetic region f3j < 1. However, the results of [5] when restricted 
to bounded degree graphs do not hold throughout the uniqueness region and hence are incomparable to ours. 
In a recent paper, Li, Lu and Yin [7] improve upon the work of [5]. They consider two-state anti- 
ferromagnetic spin systems with zero field (A = 1) on general graphs, and derive a condition under which 
an FPTAS exists for approximating the partition function. Their condition requires that (/3, 7) lies in the 
intersection of the uniqueness regions for all possible degrees d. (Note that in the (/?, 7) parameterization, 
this is a non-trivial region of the /37-plane.) However, for any fixed degree d, this region is smaller than the 
uniqueness region for d, which is the range of validity obtained in our present paper. We also point out an 
important qualitative difference between the parameterization we use (edge potential f3 and vertex field A) 
and the parameterization via two edge potentials (/3 and 7) used by Li et al: while the (/?, 7) parametrization 
is not monotonic, in the sense that uniqueness on a <i-regular tree does not imply uniqueness on a (d — 1)- 
regular tree, the (/3, A) parametrization is monotonic. In fact, as our results show, uniqueness on the 
d-regular tree implies strong spatial mixing, and hence uniqueness, on all d'-regular trees for d' < d in the 
(/3, A) parameterization. 

2 Preliminaries 

2.1 Notation 

We will mostly follow the notational conventions of [5]. Given a graph G = (V,E), a two-state spin con- 
figuration is defined as an assignment a : V i— >• {+ , — } of spins to the vertices. Weights for different 
configurations are computed in terms of the (+, +)-cdge activity j3, the (— , — ) edge activity 7 and a vertex 
activity A, and are given by 

w(a) = A m(CT) /3"+ (CT) 7 n - (CT) , (4) 

where given the configuration a, m(o~) denotes the number of vertices assigned spin — , and n+{o~) (respec- 
tively, ri_(cr)) denotes the number of edges for which both endpoints are assigned spin + (respectively, — ). 
The partition function is defined as 

Z= £ w(a). 
<?e{-i,i} v 

We remark that this representation can be easily translated to the usual description in terms of edge potentials 
and vertex field: for completeness we give the translation in Appendix [A] 

Definition 2.1 (Occupation probability). Given a vertex v in the graph G, the occupation probability 
p v is the probability that v is assigned spin + in a random configuration a sampled according to the weights 
defined in equation Q). 

2.2 The Ising model 

The Ising model corresponds to the case /3 = 7. The model is ferromagnetic when /3 > 1 and anti- 
ferromagnetic when [3 < 1 (the case j3 — 1 is trivial). The zero-field case corresponds to A = 1, the positive 
field case to A < 1 and the negative field case to A > 1. As shown in Appendix [A] on d-rcgular graphs 
the Ising model is equivalent to general two-state spin systems. Thus, in the rest of this paper, we will 



concentrate mostly on the Ising case. On non-regular graphs the equivalence still holds; however, the vertex 
activity A in the Ising model may then be different on different vertices. The adaptation of our results to 
this setting is described in Corollary |1.3| 

2.2.1 Phase transition. 

The anti-ferromagnetic Ising model exhibits a uniqueness phase transition on the d-ary tree for d > 2. In 
particular, one can define a critical activity A c (/3,d) as follows. 

Definition 2.2 (Critical activity). Consider the anti-ferromagnetic Ising model on an infinite d-ary tree 
with edge activity /3 and vertex activity X. If j3 > t=t then the Gibbs measure is unique for all values of X. 
If ft < jtt, then there exists a critical activity X c (j3,d) > 1 such that the Gibbs measure is unique if and 
only i/|logA| > log A c 08, d). 

A consequence of uniqueness^] of the Gibbs measure is weak spatial mixing, which captures a weak notion of 
decay of point to set correlations. Let p p (a, S) be the probability of occupation of the root p of an infinite 
d-ary tree when the spins of a set S of nodes are fixed according to the configuration a. Let S(p, S) denote 
the distance of p from the set S. 

Definition 2.3 (Weak spatial mixing). Given any two-state spin system, weak spatial mixing is said to 
hold if for any set S whose distance S(p, S) from the root p of the tree is finite, and any two configurations 
a i and ui, we have 

\p P (*i,S) - p p (a 2 ,S)\ < exp(-n(*(p, S))). 

Notice that weak spatial mixing does not guarantee exponential decay of correlations when the set <S* contains 
vertices which are very close to the root p, even when <j\ and o~2 differ only on vertices which are very far 
away from p. A related but, as the name suggests, stronger notion is that of strong spatial mixing, which 
captures the idea that fixing vertices near the root to the same spin should not affect the exponential decay 
of point-to-set correlations. We note that strong spatial mixing is not in general implied by weak spatial 
mixing for arbitrary spin systems; see Appendix [B] for a counterexample involving the ferromagnetic Ising 
model with appropriate parameters. 

Definition 2.4 (Strong spatial mixing). Given any two-state spin system, strong spatial mixing is said 
to hold if for any set S whose distance 5{p, S) from the root p of the tree is finite, and any two configurations 
a~\ and oi which differ only on a set T C S of vertices, we have 

\p p (a u S)-p p (a 2 ,S)\ < exp(-0(J(p,T))). 

2.2.2 Phase transition and tree recursions. 

It is well known (see, for example, [3]) that the uniqueness condition for two-state spin systems on d-ary 
trees can be written in terms of the number of fixed points of the recursion for occupation probabilities. 
Consider a subtree rooted at a vertex v in the d-ary tree, and let Uj, i = 1, 2, ...d be its children. Let p v be 
the occupation probability at vertex v and define R v — ~ Pv . One can then write the following recurrence 
for R v : 

This can easily be converted to a recurrence for occupation probabilities. Define 

A + (1-P)x 



h(x) 



l-(l-/3)x' 



As stated in the introduction, we exclude the boundary of the uniqueness region here. 



We can then write the recurrence as 

Pv=F(p Vl ,p V2 ,...,p Vd ) = - — — d — -. (5) 

We will find it useful in what follows to consider the tree recurrence with the special boundary condition in 
which all vertices at some distance I from the root are fixed to the same spin. In this case, by symmetry, 
the tree recurrence outputs the same value p v at all vertices v which are at the same distance from the root. 
Thus, the recurrence can be simplified to a one-parameter recurrence as follows: 

Pv=f(p Vl )- 



1 + Xh(p Vl ] 



Note that in the anti-ferromagnetic case, h is an increasing function, and hence F and / are decreasing in 
each of their arguments. We also note that since / is strictly decreasing in [0, 1], it has a unique fixed point. 
In terms of the recurrence function /, the condition for uniqueness can be stated as follows. 

Theorem 2.5 (|3j) For given values of ft and X, the infinite d-ary tree has a unique Gibbs measure if and 
only if the two-step recurrence function f o / has a unique fixed point. In particular, if the Gibbs measure 
is unique, and (fj, A) are not on the boundary of the uniqueness region, then the unique fixed point x* of f 
satisfies 

/V) > "I- (6) 

Remark: In [3] , it is claimed (implicitly) on the basis of numerical simulations that the condition pi is also 
sufficient for uniqueness. To be precise, the expression for the critical activity A c (/3, d) given in [3J p. 255] is 
exactly the same as that obtained by assuming that §6h is also a sufficient condition for uniqueness. While 
we believe this fact to be folklore, we have not been able to find a rigorous proof of it in the literature. With 
a slight abuse of terminology, we will henceforth refer to the set of (/?, A) for which the fixed point x* satisfies 
f'{x*) > —1 as the "uniqueness region". We will justify this terminology later (see the Remark following 
the proof of Theorem |1.1| in Section RJ by proving that condition (|6| does indeed imply uniqueness. Thus 
we will obtain a rigorous proof of the expression for the critical activity appearing in [5] . 

2.3 Messages 

Definition 2.6 (Message). A message is a continuously differentiate function (f> : [0, 1] i-> M with positive 
derivative. 

Note that a message is strictly increasing and hence invertible on its range. Moreover, the inverse function (fr 1 
is also a continuously differentiable function with positive derivative. 

Given a recurrence function / : [0, 1] H> R + , and a message </>, we denote by /^ the function ijio/o^ 1 . 
The function f^, which will play a crucial role in this paper, describes the evolution of the message 4> under 
the recurrence, in the sense that f^{4>(x)) — <fi(f(x)). We will also need the following fact. 

Fact 2.7 For any message <f), the parameters (/3, A) are in the uniqueness region if and only if f^ (p*) > —1 
at the unique fixed point p* of f^ . 

Proof. Notice that since </> is strictly increasing, and / has a unique fixed point x*, /* also has a unique 
fixed point p* — <j>(x*). Now, we notice that /^ (p*) = f'(x*), because 

/*V) = <£'(/(<rV)))/ , (<rV))<r 1 V) 

= /V). 



where in the second line we used the facts that <j>^ 1 (p*) = x* and _1 (y) = ,,, , j 1( - y. , and in the third 
line the f act t hat f(x*) = x* . Thus, ((3, A) are in the uniqueness region (as defined in the Remark following 
Theorem 



2.51 if and only if /*'(p*) = f'(x*) > -1. 



2.4 Weitz's tree reduction 

As indicated in the introduction, Weitz [T3] proved the following combinatorial reduction. 

Theorem 2.8 For any two-state spin system, strong spatial mixing on the d-ary tree implies that there 
exists a deterministic fully polynomial approximation scheme for the partition function of the spin system 
on graphs of degree at most d + 1 . 

3 Messages and contraction on the d-ary tree 

In this section, we will prove the main technical ingredient of our result, which is expressed in the following 
theorem. 

Theorem 3.1 Given d, ft and X, there exists a message <p and a constant c < 1, such that the tree recurrence 
9 — / for the quantity 4>(Pv) satisfies ||<7'||oo < c < 1 , whenever (f3,X) is in the uniqueness region for the 
d-ary tree. 

The above theorem says that in the uniqueness region, the single-parameter recurrence g = f^ for the 
message 4>(p v ) contracts at every step. (Without the message, the function / itself is not contractive.) This 
stepwise contraction is easily seen by standard arguments to imply weak spatial mixing; for completeness, we 
give a proof in Appendix |C.1| To extend the argument to strong spatial mixing, as required for Theorem |1.1[ 
we need to consider a multi-parameter (vectorized) version of the message recurrence g, since under arbitrary 
boundary conditions the occupation probabilities need not be uniform. We will show in Section HI that for 
our message <j> i n Theorem |3.1| the analysis of the vectorized version can in fact be reduced to an application 
of Theorem O 

Remark: For ease of notation, in the rest of the paper we will prove our results in terms of the uniqueness 
threshold of the d-ary tree, relating it to algorithms on graphs of degree at most d+1. As already noted, the 
uniqueness thresholds on the (d + l)-regular tree and the <i-ary tree coincide, and hence our results apply 
equally to the infinite (d + l)-regular tree. 

We begin by setting up some notation for the proof of Theorem |3.1| Notice that in the light of Fact |2.7| 



the main technical challenge is to come up with a message <f) such that the quantity 
the unique fixed point of f?. Let us fix constants 



is maximized at 



2. , n m2„„, r, VA+^-VA 



2VA 



A = d(l- P z ) + (1 - PY and D 

Define 

<f>(x) = log ( X 

Notice that D > 0, so c() is a continuously differentiable function with positive derivative on the interval [0, 1]. 
Using this message we are able to prove the following. 

Lemma 3.2 Consider the anti-ferromagnetic Ising model on a d-ary tree with edge activity j3 and vertex 
activity A. Then, defining g — f^, ip — <fi~ x , a — ip(x) and rj = /(a), we have 

a"(x) = (r, a)a'(xW(x) dftl- P)W + Aar, + A{1- a){l-r,)) 



The proof of Lemma 3.2 is somewhat technical and is given in Appendix [Dj 

Before proceeding with the technical development, we pause to give some comments on the design of our 
message. Notice that the requirement that the derivative of the function g = fv should have its maximum 
magnitude at the unique fixed point of g does not immediately lead to a solution for 0, and thus we must 
resort to some educated guesswork for the functional form of </>. Our choice is guided by the intuition 

that, by analogy with the zero field case, where it is well known that the simple message 4>{x) = log I j^ J 

is sufficient, a log ratio of probabilities shifted by an additive constant D to account for the field should 
be appropriate. The choice of D is then determined by the above requirement. An important additional 
property of our message is that, perhaps surprisingly, it does not depend upon the vertex potential A, but 
only upon the edge potential (3 and the degree d; this is reflected in the fact that the additive shift D is the 
same for both probabilities. This property will be important in extending our algorithm to the setting of 
general two-state anti- ferromagnetic spin systems in Corollary|1.3| 



Lemma 3.3 Let g = f^, with the message (f> defined above. Then \g'(x)\ is maximized at the unique positive 
fixed point of g. 



Proof. We use the notation established in Lemma 3.2 where we derived an expression for g" ' (x) in equation 
ffi\. It is easy to see that ignoring the factor (77 — a), the rest of the right hand side of equation ffi\ is 
negative: this is because g is a decreasing function, while ip, being the inverse of the increasing function cf), 
is increasing. Also, we have < f3 < 1 (in the anti-ferromagnetic case) and < a, rj < 1 (since they are 
probabilities), so that the fractions appearing on the right hand side are positive. 

Let x* be the unique fixed point of the strictly decreasing function g. From the above discussion, it 
follows that the sign of g"(x) is the opposite of the sign of n — a = f(ip(x)) — ip(x). Notice that r\ — a is 
strictly positive for x < x* and strictly negative for x > x* . This implies that g'(x) is strictly decreasing for 
x < x* and strictly increasing for x > x* . Since g is strictly decreasing this shows that the magnitude of g' 
is maximized at a;*. I 



Combining Lemma 3.3 with Fact|2.7| we immediately get Theorem 3.1 Lemma 3.3 further implies that the 



constant c in the Theorem is |</(:r*)|, where x* is the unique fixed point of g. 

4 Strong spatial mixing on the d-ary tree 

In this final section we use the message defined in the previous section to prove our main result, Theorem |1.1| 
Along with Weitz's reduction stated in Theorem |2.8[ this will immediately imply Corollary |1. 2 1 the FPTAS 
for the anti-ferromagnetic Ising model with arbitrary fields. To derive Corollary |1.3| for general two-state 
spin systems, we will need the translation described in Appendix [A] Both these latter proofs appear at the 
end of this section. 

Recall that Theorem |1.1| asserts that weak spatial mixing implies strong spatial mixing. We already 
showed in Theorem |3.1| that in the uniqueness region, there is uniform contraction in the tree recurrence 
with uniform boundary conditions. However, in order to prove strong spatial mixing, we will need to handle 
non-uniform boundary conditions as well, in which case the one-parameter recurrence g = /* is no longer 
sufficient. We therefore consider the multi-parameter vectorized version G of the function g. For x € M. d , 
G(x) is defined as 

G{xx,X2,...,x d ) = <t> (- — — d — — — - 1 , (8) 

where, as before, ip — <f>~ 1 ■ We claim that strong spatial mixing on the d-ary tree is implied whenever the 
function G satisfies the following condition; a proof of this implication can be found in Appendix |C. 2 1 

Definition 4.1 (Contractive spatial mixing). Given the parameters j3 and A for the anti-ferromagnetic 
Ising model on a d-regular tree, contractive spatial mixing holds if there exists a constant c < 1 such that 

\G{x)-G{y)\<c\\S-y\\oo, 
10 



for the vectorized version G of g defined as above with respect to the message <fi. 

To establish this condition, we will rely on the following lemma. 

Lemma 4.2 Let r) = ip(G(x)). Let fj be the unique solution of il>(g(ff)) — 1). Then ||V(G(x))||i < 
\\V(G(n,fj,...,fj))h = W(v)\- 

Proof. Set on = ip{xi) for i = 1, 2, . . . d. We then have 

??= l + AIltiM«i) = l + A/i(^)) d ' (9) 

Recalling the definitions of the quantities A and D given just before Lemma |3~2J we can now write ||V (G(x)) ||i 
as 

\\*7fnf#\\\\ dn(l - n)(l - (3 2 ) ( 2 ^ a. t (l-«. t ) \ 

llv(o(*))lli- p + M1 . v) ( 1 + ( 1 -^)g^ + (1 _ /3)2ai(1 _ ai )J- ( 1Q ) 

For notational convenience, we define the function J(x) = g,( 1 _g\2 X i 1 _ \ ■ Note that maximizing the sum in 






10) under the constraint rt9t is the same as maximizing Xh=i J( a i) under the constraint that Y\i=i h(°ti) = 
j^ 1 . Since /i is positive and invertiblc, it is therefore sufficient to show that the function K(x) = J(h^ 1 (e x )) 

is concave in order to show that all aj's are equal at a maximum. We now show this by direct computation. 

After differentiating twice and simplifying, we have 

K"(x) = - e ~ X{1 + e2X) P <0 



This shows that K is concave. By the discussion above, it follows that the sum in equation ( 10 ) is maximized 
when all ctj's are equal. In conjunction with the condition that r\ = - — d l — r, this shows that 

||V(G(f))|| 1 <||V(G(r7,,7,...,^))|| 1 . 
Note that for any x, G(x, x, ...,x) = g(x), and therefore ||V (G(fj,fj, . . . ,fj))\\i ~ \g'(fj)\. I 



Using Lemma 3.3 and the above lemma, we are now ready to prove our main technical result, Theorem |1.1[ 



which says that weak spatial mixing implies strong spatial mixing for general anti-ferromagnetic Ising models. 

Proof of Theorem \l.l\ Consider a setting of parameters j3 and A such that the d-ary tree has weak spatial 
mixing. Let x* be the unique fixed point of the function g. We will use only the property that the fixed 
point satisfies the condition (pi of Theorem 2.5 By Theorem 3.1 we have |l<?'||oo = c < 1. By Lemma 4.2 



this implies that for all x in the domain of the function G defined in equation pi), ||VG(a;)||i < c. Using the 
mean value theorem followed by Holder's inequality, we then have 

\G(x) - G(y)\ < c\\S - yWoo, 

for all vectors x and y in the domain of G, and thus contractive spatial mixing holds. As observed above, 
this implies strong spatial mixing. I 

Remark: We can now justify our use of the term "uniqueness region" as described in the Remark following 
Theorem |2.5| Notice that in the proof of Theorem |1 . 1| above, we used only the fact that weak spatial mixing 
implies that (/3, A) is in the "uniqueness region" as defined in the aforementioned Remark. Thus, we see 
that whenever (/3, A) is in the uniqueness region, we have strong spatial mixing, and hence, in particular, 
uniqueness. As stated earlier, this provides a rigorous proof of the claim in [3] that the interior of the 
uniqueness region is equivalent to the condition ([6]). 



we can now 



Combining the above theorem with the general reduction of Weitz [TJ] stated in Theorem 2.8 
prove Corollary [L2j which asserts the existence of an FPTAS for general anti- ferromagnetic Ising models on 
bounded-degree graphs up to the uniqueness threshold. 
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Proof of Corollary \1.£\ As observed earlier, in order to obtain an FPTAS for the partition function of the 
associated spin system, it is sufficient to give an FPTAS for approximating the occupation probability p p of 
a vertex p, under an arbitrary fixing of spin values for an arbitrary subset of vertices. Given a vertex p in a 
graph G of maximum degree (d+ 1), we start by constructing Weitz's self-avoiding walk (SAW) tree rooted 
at p. For non-leaf vertices (apart from p) in this tree which do not have d children, we can create dummy 
children (so as to make the arity of the vertex d) all of which independently have occupation probabilities 
of 1/2. It is easy to see that this does not change the output of the tree recurrence (equation ^) at any 
vertex of the tree. As we saw in the proof of Theorem |1.1| we have strong spatial mixing on this SAW 
tree whenever (/3, A) are in the uniqueness region of the d-ary tree. The corollary now follows using Weitz's 



reduction (Theorem 2. 



Finally, we will see how to use Lemmas |3.2| and |4.2| to prove Corollary |1.3| which extends the FPTAS to 
general two-state anti-ferromagnetic spin systems. 

Proof of Corollary \1.3[ Given a two-state spin system with parameters /?, 7 and A on a graph G of degree 
at most d+ 1, we can use the translation given in Appendix [X] to come up with an equivalent Ising model 
with edge potential /3' = \[Wl and vertex-dependent potentials A„ = \(y/"//fl) d "" . Now, as before, in order to 
estimate the occupation probability p p for a given vertex p, we construct Weitz's self-avoiding walk (SAW) 
tree rooted at p, and complete the degree of any non-leaf vertex (apart from p) in the tree which does not 
have d children by attaching dummy children which are fixed to have occupation probability | . We now use 
the message <f> constructed above for <i-ary trees for the parameter /3'. By the hypotheses of the corollary, the 
parameters (j3', X u ) at each vertex u of the SAW tree are in the uniqueness region of the d-ary tree. Since 
the message <f> does not depend upon \ u , Theorem |3.1| and Lemma |4.2| apply at each vertex u of the tree. 
Thus, as in the proof of Theorem |1.1| we get contractive spatial mixing and, hence, strong spatial mixing 
on the SAW tree. Employing Weitz's reduction (Theorem |2.8[ ), we have the first part of the corollary. 

The claim that the class Q in the corollary includes id + l)-regular graphs when f3, 7 and A are in the 
uniqueness region of <i-ary tree follows by noticing that in this case the parameters A' = X v obtained by the 
translation are the same at each vertex v, and that j3' and A' are in the uniqueness region of the d-ary tree 
by the hypotheses of the corollary. Thus, we can complete the proof for this case in the same manner as in 
the proof of Corollary |1.2[ I 
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Appendix 

A Translation between various descriptions of the Ising model 

General two-state spin systems are usually described in terms of (symmetric) energy functions Q(+ , + ), 
Q(+ , — ) = Q(— , + ) and Q{— , — ), and an odd vertex field h(+) — —h(—) — h. For a graph G — (V, E), 
the partition function of the system is then Z 2 — J^ w 2 (c) , where the sum is over all states of the system 
a : V —t {+, —} and w<2{a) is defined as 

u> 2 (cr) = exp I - Y^ Q(< J ( u )'< T ( v ))-^2 h ( (J ( v )) 

\ {u,v}eE v£V 

This is in fact equivalent to our formulation of the system given in Section |1.1| (see equations (IT]) and (pi)) 
To see this, define 

/3 = cxp(-Q(+,+) + Q(+,-)); 
7 = cxp(-Q(-,-) + Q(+,-)); 
A = exp(2/i), 



which yields 

for all a, and, therefore, 



w(a) = w 2 {a) exp (Q(+, -)\E\ + h\V\) 



Z = Z 2 c^{Q{+,-)\E\ + h\V\). 



We call the above spin systems soft constraint systems if B, 7 and A are non-zero, or equivalently, if the energy 
functions and field are finite for all spin values. As we shall now see, every such soft constraint system can be 

13 



represented in terms of the Ising model (this translation can also be found, e.g., in [S]). Consider a general 
two-state spin system with parameters /?, 7 > and A. Then the equivalent Ising model has edge activity 



and a degree-dependent vertex activity given by 

A., = 




where d v denotes the degree of vertex v. Now, denote the weight of a configuration a in the Ising model just 
defined by w*(a) and its partition function by Z*. Then one calculates straightforwardly that 

w(a) — w*(a) 
and hence 





Thus we have translated the original spin system with parameters (/3,7, A) into an Ising model with 
locally changing field. Note that on regular graphs the resulting field is in fact constant at all vertices. 
Furthermore, the Ising model is anti- ferromagnetic if and only if /?7 < 1. This justifies our use of the term 
"anti-ferromagnetic" for general spin systems based on the value of fij. We also observe that in the special 
case of d-regular trees, this implies that weak (strong) spatial mixing in the original spin system (/3, 7, A) 
is equivalent to weak (strong) spatial mixing in the Ising model given by the translation. A little thought 
shows that since all vertices except the root have the same degree in the d-ary tree, the last observation 
holds also for d-ary trees. 

B Weak spatial mixing does not imply strong spatial mixing for 
the ferromagnetic Ising model 

We construct a counterexample as follows: given a degree d > 3, consider the infinite rooted d-ary tree, with 
the fixed boundary condition where each vertex in the tree has one of its children fixed to +. Notice that 
if the original parameters are /3 > 1 and A, the effect of this fixed boundary condition can be simulated 
by changing the vertex field to -5 . Therefore, strong spatial mixing on this subgraph of the d-ary tree with 

parameters (/?, A) holds only if weak spatial mixing holds on the id — l)-ary tree with parameters (f3, h). It 
is therefore sufficient to choose j3 and A satisfying both the conditions 



log A > log A c (/3,d), and (11) 

0<logA-log/3<logA c (/3,d-l) (12) 

in order to construct a counterexample. To see that such a choice of parameters is possible, we consider the 
exact form of X c (f3, d). Translating the results in [3} p. 250] to our notation using Appendix [A[ we have 

log\ c (p,d) = (d-l)logp-P(d) + Q(p,d), 

where 

P(d)=dlogd-(d-l)log(d-l), and 

lim Q(j3,d) — 0, for any fixed d. 
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We note that P(d) is an increasing function of d. Thus, the required conditions (11) and (12 1 become 



logA>log/3, (13) 

logA> (d-l)log/3-P(d) + QCM), and (14) 

logA< (d-l)log(3-P{d-l) + Q((3,d). (15) 



For a fixed d, Q{f3,d) is 0/3(1), and hence inequality (14) implies inequality (13) for /3 large enough and d > 3 



Again, since Q{(3,d) is 0/3(1), and P{d) is an increasing function, it is possible to find A satisfying both the 



inequalities (TT4J) and (15) for /3 large enough. Thus, for any d > 3, we can find (J3, A) with /3 > 1 such that 



weak spatial mixing for the d-ary tree does not imply strong spatial mixing. 

C Contraction and spatial mixing 

C.l Contraction and weak spatial mixing 

We show in this section that a stepwise contraction in the recurrence for 4>{p v ) implies weak spatial mixing. 
As before, we denote /^ by g, and assume that for any x, y in the range of <fi, it holds that 

\g(x) - g(y)\ < c\x - y\, (16) 

for some c < 1. To show that this implies weak spatial mixing, we consider boundary conditions u\ and 02 
on a set S whose distance from the root p is I. Using the monotonicity of the tree recurrence F (defined in 
equation (|5|)) in all its arguments, it can be verified that \(j>(p p (ai, S)) — (j)(p p (a2, S))\ is maximized when S 
is the set of all leaves at distance I from p and o\ assigns all vertices in S to + and <Ji assigns all vertices 
in S to — . With this definition of a± and o"2, we notice that the tree recurrence for <fi(p v ) outputs the same 
value for all vertices v at the same distance from the root. For a vertex at distance I — i from the root p, we 



denote by qi j the quantity 4>(p v {cr j , S)) , for j £ {1,2}. Now, using condition (16), we have 

ki+1,1 - Qi+iM = \g(qi,i) - g(qi,2)\ 

< c\qi t i - q t , 2 1 ■ 

Since both 4> and (j>~ 1 are continuously difierentiable functions defined over compact sets, they are Lip- 
schitz continuous, say with parameters L\ and L2 respectively. We therefore have weak spatial mixing, 
since 

\p p ((Ji,S) -p p (a 2 ,S)\ = |0 _1 (g ; ,i) -0 _1 (g;, 2 )| 

< L 2 c l \q ,i -qaal 

< L x L 2 c l . 

C.2 Contraction and strong spatial mixing 



In this section, we show that contractive spatial mixing, as defined in Definition 4.1 implies strong spatial 
mixing. We again consider boundary conditions o\ and a 2 on a set S which differ only on a subset T which 
is at distance I from the root p. Again, since both <f> and </> -1 are continuously diffcrentiable functions defined 
over compact sets, they are Lipschitz continuous, say with parameters L\ and L 2 respectively. We define 
the quantity q L as 

qi= max \(f>(p v (ai, S)) - <fi(p v (a 2 , S))\ . 

v:S(p,v)=l — i 

Notice that go < 10(1) — 0(0)1 ^ -^l- Also, since G is the tree recurrence for <fi(p v ), contractive spatial 
mixing for G implies qi+\ < cqi for c < 1. Thus, we get strong spatial mixing since 



\p p (ai,S) ~p p (a 2l S)\ < L 2 qi < L 2 c l q < LiL 



>2C>. 



15 



D Proof of Lemma 13.2 



In this section, we prove Lemma |3.2| The proof involves a few somewhat lengthy derivative computations, 
which we isolate in the following lemma. 



(17) 
(18) 
(19) 
(20) 
(21) 



Lemma D.l With the notation used in Lemma \3.S\ above, we have 

cf>"(x) _ A(2x-1) 

4>'{x) ~ fi + Ax{l-x)' 
h'{x) _ 1-/3 2 



h(x) /3 + (1 - /3) 2 x(l - x)' 1 
h"{x) _ 2(1 - 0) 
h'(x) ~ l-(l-/3)x' 



f(x) = -df(x)(l-f(x)) 



h'jx) 
h(x)'' 



f"(x) = f(x)(l-2f(x)) h"(x) h'(x) 
fix) f(x)(l-f(x)) h'(x) h(x)- 



Proof (sketch). Most of these identities are easily verified by direct computation. In proving equation (17 1, 
one needs to keep in mind the definition of the constant D. I 



Proof of Lemma \3.ty To ease notation, we will suppress the dependence of the quantities 77 and a on x. 
Using the chain rule, we have 



9 ( x ) = TTTTT/ (")• 



«/)'(«)• 



Here, we used the fact that since tp = (j> l , ifj'(x) 



After taking the logarithm, and noticing that 



the right hand side is more easily expressed as a function of a rather than of x, one can write the second 
derivative of g as 

•4>'{x) g'{x) <j)'{rj) da 4>'{a) /'(a) ' l ' 



We now consider each of the terms involved above. Recalling that r\ = /(a), and using equations (20) 



and (21 1 to expand the first and last terms in equation (22) above, we get 

1 fix) 



ift'ix) g'ix) 



= T 1 -T 2l 



(23) 



where T\ and T2 arc defined as 



T, 



h"ja) h'ja) (j>"{a) 
h'(a) h(a) <P{a) 



, and 



T 2 ^d 



h'ja) 
h(a) 



<f>"(v) 
<f>'(v) 



77(1 — 77) + 1 - 2?/ 



Notice that all terms containing n are isolated in T 2 . We now consider each of the terms separately. For Ti, 

we have 



h"ia) h'ia) _ 2(1 - 0) 



1-P 2 



h'ia) hia) 1 - (1 - 0)a /3 + (1 - ^) 2 a(l - a) 
(l-/? 2 )(2a-l) 
j3 + (1 - /3) 2 a(l - a) ' 



16 



Here, we used equations (Tl9|) and (18 1 in the first line. Now using equation (17), we have 



Ti = 



(2a - 1) ((1 - jif \p + Aa(l - a)} - A [/? + (1 - pfajl - a)] ) 
(ft + (1 - /3) 2 a(l - a)) (/? + Aa(l - a)) 

/3(2a-l)((l-/3) 2 -A) 
09 + (1 - /3) 2 a(l - a)) 08 + Aa(l - a)) 

-dp(2a-l) /i'(a) 



(/3 + A*(l - a)) /i(a) 



Here, we use A = d(l — /3 2 ) + (1 — /3) 2 , followed by equation (18) in the last line. 
We now consider T%. Again using equation (17), we have 



T 9 



h'{a) 



^(277-1)77(1-77) 



h(a) [ P + An{l -77) 
-d/3(277- 1) ti(a) 
P + Ar)(l - rj) h{a) ' 



- (2r? - 1) 



Notice that modulo the 



h'(a) 



factor, T\ and T 2 have the same functional form as functions of a and 77 



respectively. In fact, the message <f> is designed so as to make this possible. We can now substitute these 

277 — I 2a - 1 



values into equation ( 23 ) to get 



g"(x) = d(3g'(xW(x 



ti{a) 
h(a) 



dpg'{xW{x) 

(ri-a)g'(x)ip'(x) 



/3 + Ar)(l-r)) P + Aa(l - a) _ 
ti(a) (77 - q)(2/3 + Ajar] + (1 - q)(l - 77))) 
h(a) {0 + Aa(l - a)) {fl + Aq(l - rj)) 

d/3(l - /3 2 )(2/3 + Aar, + A{1 - a)(l - 77)) 



(p + o(l - a)(l - /3) 2 )(/3 + Ar?(l - r;))(/3 + Aa(l - a)) ' 



where in the last step we used equation (18) 
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